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Theoretical difficulties in reconciling the measured rates for ordinary and radiative muon capture are 
discussed, based on heavy-baryon chiral perturbation theory. We also examine ambiguity in our analysis 
due to the formation of pfip molecules in the liquid hydrogen target. 

1 Introduction 

o 

Ordinary and radiative muon captures (OMC and RMC) on a proton, fi~p — » nv^ and fi~p — ► nv^j, 
are fundamental weak-interaction processes in nuclear physics and constitute primary sources of 
information on gp, the induced pseudoscalar coupling constant of the weak nucleon current [lj. 
The most accurate existing measurements of the OMC and RMC rates have been carried out using 
a liquid hydrogen target. The experimental OMC rate in liquid hydrogen obtained by Bardin et 
at. 01 is 



>' = 460 ± 20 [s- 1 ] (OMC). (1) 

As for RMC, Jonkmans et al. 0] measured the absolute photon spectrum for > 60 MeV and 
deduced therefrom the partial RMC branching ratio, i? 7 , which is the number of RMC events (per 
stopped muon) producing a photon with E y > 60 MeV. The measured value of i? 7 is f3L |j 

R^p = (2.10 ± 0.22) x 10~ 8 (RMC). (2) 

Surprisingly, the value of gp deduced in |3| from the RMC data is ~1.5 times larger than the 
PCAC prediction, gp CAC ||. By contrast, the value of gp deduced in Q from the OMC data is in 
good agreement with the PCAC prediction. Heavy-baryon chiral perturbation theory (HBChPT), a 
low-energy effective theory of QCD, allows us to go beyond the PCAC approach, but the results of 
detailed HBChPT calculations up to next-to-next-to-leading order ||] essentially agree with those 
obtained in the PCAC approach. Thus the theoretical framework for estimating gp appears to be 
robust. What then can be the origin of the apparent conflict between the gp values determined 
from OMC and RMC ? In this talk we wish to address a number of issues relevant to this question. 

The OMC process is described by the standard electroweak theory. The lepton current, j l ^P ton ^ 
interact with the very heavy W boson (mw — 80 GeV) which propagates and interacts with 
the hadronic current, jj^ adron . Since the maximal momentum transfer between the two currents, 
q max ~ <C mw, we can for all practical purposes replace the W boson propagator with a 
constant. Then the interaction between the lepton and hadron currents reduces to: 

Gf -lepton jhadron /o\ 
~ -/= Ja J a ■ ( 3 ) 

where Gf is the effective Fermi constant. Since the leptons are considered point particles, the lepton 
current is given by the incoming muon and outgoing neutrino spinors: j l ^P ton = u^a^l — 75)^. 
The hadronic current however is less well known due to the structure of the nucleons. We write 
jhadron _ _ ^ ^ where the vector current, Vp, and the axial current, Ap, can be written as 



(based on symmetry considerations and neglecting second-class currents): 



Vp = u r , 
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where u n and u p are the outgoing neutron and incoming proton spinors; and q is the momentum 
transferred to the nucleon. The four form factors, gv{q 2 ): 9m{<1 2 ), 9a{<1 2 ) and gp{q 2 ), in Eqs.@ 
and (|) are, at low momentum transfers relevant to the OMC and RMC reactions, given by known 
parameters: 

9v(q 2 ) = l + \<r 2 >q 2 ^ ; gu{q 2 ) = k p - K n H ; (6) 

b 

9A(q 2 ) = 9a(i + \ <r\>q 2 + --^j ; (7) 

9p{q ) = — ^ n odAm^iriN <r A > H . (8) 

m£ — 3 

Here < r 2 > ~ 0.585 fm 2 is the square of the isovector nucleon radius, < r\ > ~ 0.42 fm 2 is the 
axial radius squared, and n p and K n are the proton and neutron anomalous magnetic moments, 
respectively. In the induced pseudoscalar form factor, gp(q 2 ), the form of the dominant "pion-pole" 
term is dictated by the PCAC hypothesis. 

In OMC the momentum transfer is q 2 =(m fl — E v ) 2 — E 2 =— 0.88m 2 and, traditionally, the 
pseudoscalar coupling "constant", gp, is defined as gp = gp(q 2 = —0.88m 2 ), and thus the PCAC 
value is gp CAC = 6.77^(0) We remark that the contribution of the gp(q 2 ) term to OMC is 
suppressed due to the fact that q 2 =— 0.88m 2 is far from the pion pole. The rationale for using 
RMC to determine gp is that the three-body kinematics in the final state allows q 2 to approach the 
pion pole, enhancing thereby the contribution of the pseudoscalar term. With the photon energy 
denoted by E 1 , we have q 2 ~ (2m fJi E ry — m 2 ), which is positive for sufficiently large values of E 1 , 
and which can reach the maximal value ~ m 2 . In ||, |j, the RMC process is measured for E^ > 
60 MeV. Meanwhile, a great challenge in observing RMC is its very low branching ratio (~ 10~ 3 
compared to OMC), and it was quite a feat that the TRIUMF group succeeded in measuring FLy. It 
should also be mentioned that the absolute photon spectrum (for > 60 MeV) measured in 
gives information on the q 2 dependence of gp(q 2 )- The afore-mentioned surprising conclusion that, 
to fit the TRIUMF RMC data, one needs to adopt gp ~ 1.5 gp CAC has triggered re-examination 
of the foundation of the theoretical framework used for the analysis. In particular, it motivated 
detailed systematic calculations based on HBChPT for RMC as well as OMC §, |, |l|, 0, 

2 Chiral perturbation theory and the atomic capture rates 

Chiral perturbation theory (ChPT) is an effective field theory of QCD tailored for describing low- 
energy hadronic interactions. Since the physical degrees of freedom here are hadrons, not quarks 
or gluons, we "integrate out" the quark and gluon fields to obtain an effective lagrangian, C e ff, 
pertinent to the hadronic degrees of freedom. C e ff should inherit all the symmetry properties of 
QCD (and the patterns of their breaking, if any), including chiral symmetry. For our purposes, it 
is sufficient to consider only two flavors (u and d); correspondingly, C e ff contains only the nucleon 
and pion fields. If chiral symmetry is an exact symmetry of £ e ff: then the right- and left-handed 
nucleon fields decouple, and the left- and right-handed Noether currents are separately conserved: 
d a jf e j t = and d a j" ight = 0. The QCD vacuum state, however, does not respect chiral symmetry, 
i.e. chiral symmetry is spontaneously broken. Then, according to the Goldstone theorem, there 
appear massless pseudoscalar Goldstone fields, which can be identified with (massless) pions. In 
reality, Cqcd contains a tiny quark mass term (m quar k —10 MeV) that violates chiral symmetry 
and that gives rise to a finite pion mass. A consequence of the non-zero mass is that the axial current 



is no longer conserved. Since m quar k <C h-chirai — 1 GeV, this explicit chiral symmetry breaking 
can be treated as a small perturbation, and ChPT consists in expanding C e ff an d transition 
amplitudes in terms of Q/ ^-chiral an d iTipion/ ^-chiral, where Q is a typical scale of external momenta: 
Q ~ l^pionl ~ \Pnucieon\- A version of ChPT that is particularly useful for our present purposes is 
HBChPT, wherein the nucleon is treated as a heavy particle for which a Foldy-Wouthuysen-like 
non-relativistic expansion can be used. In HBChPT we have an expansion 
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]T£ /3 = £i+£2 + 



0) 



Here the "chiral index" f3 is given by (3 = d + n/2 — 1, where n is the number of nucleon fields 
involved in a given vertex and d the number of derivatives or powers of involved. For example, 
the lowest chiral order lagrangian is given as 



d = N 

while the next order Lagrangian is given by 

C 2 = N 
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It is to be noted that in HBChPT the Schroedinger operator V 2 /(2mjv) for the nucleon kinetic 
energy is treated as a "recoil" correction to C\. In muon capture (both OMC and RMC), Q ~ 
m^=105.7 MeV or Q/A x ~ 0.1 and hence the chiral expansion is expected to converge rapidly. The 
explicit calculations in HBChPT || |l^] corroborates this expectation. 

For OMC, leading-order (LO) contributions in chiral counting come from two tree diagrams of 
order (Q / A c hirai) C '• Next-to-leading order (NLO) contributions are again given by two tree diagrams 
which however are of order {Q / K c hi ra i) 1 . These NLO diagrams arising from £2 define the nucleon- 
weak current (or pion) vertices which include the nucleon recoil correction of order ~ mj^ 1 . To 
next-to-next-to-leading order (NNLO), (Q / A-chirai) 2 , we have both tree and one-pion-loop diagrams. 
The tree diagrams here involve vertices coming from £3, which contains three low energy constants 
(LEC). Two of them can be determined from < r 2 > and < r\ >, while the third LEC can 
be constrained by the Goldberger-Treiman discrepancy. The one-pion-loop diagrams effectively 
introduce a form factor at the nucleon-pion vertex. It should be stressed that there are no free 
parameters in this ChPT calculation; all the LEC's are given by gA, k p — n n , f n ~ 93 MeV, flViv, 
<r 2 >, <r A >, and the Goldberger-Treiman discrepancy defined by m]ygA = /vr57rAf(l — A^-). 
Bernard et al. |11] have shown the very rapid convergence of these ChPT calculations for the spin- 
averaged OMC rate: A = (247 — 62 — 4 + • • •) s ~ 1 ) where the the first, second and third terms 
correspond to the LO, NLO and NNLO contributions, respectively. We can see that the NLO term 
is ~25% of the LO term and that the NNLO term is ~1.6% of the LO term. Meanwhile, the OMC 
rate from the fi-p singlet atomic state exhibits a more subtle convergence behavior, indicating that 
a systematic calculation based on well-defined chiral expansion is indeed needed to get accurate 
predictions: 
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The second term in this expression is the nucleon recoil correction. Of the two terms in the square 
brackets representing the NNLO corrections, the first represents the m^ 2 correction term while the 
second term originates from the g 2 -dependence in the hadronic form factors given by the ChPT 
terms of one-loop order. It is noteworthy that these two terms cancel each other to a significant 
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Table 1: Comparison of calculated atomic OMC rates. A s (A t ) is the atomic capture rate (in sec 1 ) 
calculated for the initial singlet (triplet) hyperfine state including terms up to NLO or NNLO, as 



indicated. The entries for the columns labeled "BHM" and "AMK" are taken from [11] and [10], 
respectively. The "AMK" results have been obtained with the use of g& = 1.267 and g n pj = 13.4. 
Apart from small chiral corrections, the numerical results of the classic works by Primakoff and 
Opat would be close to those of AMK(NNLO), if the updated value of gA = 1.267 is adopted; 
Primakoff and Opat used gA = 1-24 and gA = 1.22, respectively 



degree, a feature that cannot be studied reliably without a systematic expansion scheme such as 
ChPT. We compare in Table |l] the OMC rates obtained in two independent ChPT calculations, 



BHM |ll| and AMK @ . The variance in the results of BHM and AMK are ascribable to the slight 
differences in the approaches used. 

The RMC atomic capture rates have been calculated with the same £ e ff- The use of £ e ff 
ensures that the photon coupling in our RMC calculation automatically satisfies gauge invariance. 
The leading order (LO) diagram representing the emission of a photon by a muon gives a dominant 
contribution. In LO we also have Feynman diagrams of the following types: (i) a photon couples 
to an intermediate pion propagator; (ii) a photon couples to a pion-nucleon vertex; (iii) a photon 
couples to a H^~-pion vertex. In the Coulomb gauge, there are five LO Feynman diagrams, ten 
NLO diagrams, and more than twenty NNLO diagrams including pion-loop diagrams ||. 

It has been pointed out [|l4|] that, with the use of the atomic RMC rates calculated model- 
independently with the use of HBChPT, it is extremely difficult to reproduce Rl xp obtained in the 
TRIUMF experiment. We describe below some salient features of this difficulty. 

3 Muonic states in liquid hydrogen 

To make a comparison between theory and experiment, one needs to relate the theoretically cal- 
culated atomic OMC and RMC rates to the capture rates measured in a liquid H2 target, An q 
and i? 7 , respectively. It is also important to know the temporal behavior of each of the various 
//-capture components (capture from the atomic states and capture from p-fi-p molecular states). 
Fig. H schematically depicts various competing atomic and molecular processes in liquid H2- A 
muon stopped in liquid hydrogen quickly forms a muonic atom (/i-p) in the lowest Bohr state. The 
atomic hyper fine-triplet state (S=l) decays extremely rapidly to the singlet state (S=0), with a 
transition rate A10 — 1-7 x 10 10 s . In the liquid hydrogen target a muonic atom and a hydrogen 
molecule collide with each other to form a p-fi-p molecule with the molecule predominantly in its 
ortho state. The transition rate from the atomic singlet state to the ortho p-\x-p molecular state, 
Xppfj,, has an averaged value A PPM ~ 2.5 x 10 6 s , which is comparable to the muon decay rate, 
Ao = 0.455 x 10 6 s _1 . The ortho p-\i,-p state further decays to the para p-^-p molecular state with 
a rate A op ~ (4 - 7) x 10 4 s -1 . 

We now discuss very briefly the time structure relevant to the OMC experiment (2j . We denote 
by N s (t), N om (t) and N pm (t) the numbers of muons at time t in the atomic singlet, ortho-molecular, 
and para-molecular states, respectively. They satisfy coupled differential equations (the kinetic 
equations), see Eq. (54a) in Ref. [fl3f| . For illustration purposes, let us consider a case in which there 
is one muon in the atomic singlet state at t = 0; i.e., N s (0) = 1 and iV om (0) = N pm (0) = 0. N s (t), 
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Figure 1: Atomic and molecular states rel- 
evant to muon capture in liquid hydrogen; 
Xppfj, is the transition rate from the atomic 
singlet state to the ortho p-fi-p molecular 
state, and X op is that from the ortho to para 
molecular state. 
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Figure 2: Number of muons at time t in each 
state in liquid hydrogen. 



N om (t) and N pm (t) corresponding to this case are plotted in Fig. ||. It turns out to be crucially 
important to take proper account of the t-dependence of these populations in analyzing the OMC 
data in ||, since data taking in Ref. || starts at t ^ 0. 

In the OMC experiment (see Fig. 4 in Ref. j|]), \i~ beams arrive at the target in a (on the 
average) 3 /is-long burst with repetition rate 3000 Hz. The data collection typically starts 1 fis 
after the end of the 3 iis-long beam burst, and the measurement lasts until 306 /is after the end of 
the beam burst. To proceed, we assume that the quoted average time intervals represent the actual 
values (ignoring fluctuations). Then, provided all the muons arrive at the same time, we can choose 
with no ambiguity that arrival time as the origin of time (t = 0) and let t = ti, the starting time for 
data collection, refer to that origin. However, the finite duration (3 /xs) of the beam burst causes 
uncertainty in the value of t = fy; ti can be anywhere between 1.0 /is and 4.0/is. We choose here to 
average over the muon burst duration time for deducing An q , see Ref. [|14|] for details. 

For the RMC experiment || the muons essentially arrive one by one and data taking begins 
at ti = 365 ns. We therefore can neglect the beam burst duration time in calculating i? 7 , see 



Ref. g4|. 



4 Discussion 

The value of An g obtained with the use of the atomic OMC rates calculated in HBChPT up to 
NNLO is Aug ~ 459 s -1 |10, |j~4| . This value is in good agreement with A^ p . By contrast, it is not 



possible to reproduce R^ xp in Eq.(|l|) in the existing theoretical framework and with the use of the 
standard set of input parameters, see Ref. [14] for a more detailed discussion. 

The important question is: Have we exhausted all possibile ways for reconciling the measured 
OMC and RMC rates? One thing worth studying as a speculative possibility [|14|] is the sensitivity 
of Aug and i? 7 to a so-far neglected possible change in the value of the molecular mixing parameter 
£. As first discussed by Weinberg |I5| , £ parametrizes a possible mixing of the spin-3/2 and spin-1/2 
states in the p-jJL-p molecule, and this mixing changes the molecular capture rate to 



Km = £A om (l/2) + (1 - £)A om (3/2), 



(13) 



where A om (l/2) = 
0.25A t ) and 2 7o 



A om and A om (3/2) = 27 A 4 for both OMC and RMC; A om = 2 7o (0.75A s + 
= 1.009 [jnj. Although the existing theoretical estimate suggests £ ~ 1 |i"3f , 



we treat £, as we did in Ref. JT^J , as a parameter to fit the data. Our study [14] demonstrates 
that, even with this extra adjustable parameter, it is impossible to fit the OMC and the RMC 



data simultaneously. In Ref. 1 14] we have also found that both the OMC and RMC capture rates 



are sensitive to the value of \ op . Obviously, the results of a more precise measurement of X op at 



TRIUMF |J will be very important for both OMC and RMC. Note that since HBChPT constrains 
the value of gp with high accuracy, there is not much room for adjusting the value of gp. 

Our findings reported in |l4j and briefly summarized in this talk are mostly the reconfirmation 
of the conclusions stated in one way or another in the literature, but we hope that the coherent 
treatment of OMC and RMC in liquid hydrogen as described in Ref. Jl4|] would be useful. Although 
we have presented examples of simulation of the experimental conditions, they are only meant to 
serve illustrative purposes. Definitive analyses can be done only by the people who carried out the 
relevant experiments. Finally, we remark that a precise measurement of the OMC rate in hydrogen 
gas is planned at PSI[16|. This experiment will allow us to avoid the molecular complexity discussed 



above and directly test the HBChPT prediction [1C, 11 1 



Acknowledgments: We are grateful to Peter Kammel, Tim Gorringe and Harold Fearing for useful com- 
ments. This work is supported in part by NSF grant PHY-9900756. 



References 

[1] See, e.g., M. Morita, Beta decay and muon capture, W. A. Benjamin, Inc. (1973). 

[2] G. Bardin et al, Phys. Lett. B 104 (1981) 320; Nucl. Phys. A 352 (1981) 365. 

[3] G. Jonkmans et al, Phys. Rev. Lett. 77 (1996) 4512. 

[4] D. H. Wright et al, Phys. Rev. C 57 (1998) 373. 

[5] See, e.g., S.L. Adler and Y. Dothan, Phys.Rev. 151,1267 (1966) 

[6] V. Bernard, N. Kaiser and U.-G. Meifiner, Phys. Rev. D50, 6899 (1994); V. Bernard, H.W. Fearing, 

T.R. Hemmert and U.-G. Meifiner, Nucl. Phys. A 635, 121 (1998); ibid A 642, 563 (1998) 
[7] H. W. Fearing, Phys. Rev. C 21 (1980) 1951. 

[8] T. Meissner, F. Myhrer and K. Kubodera, Phys. Lett. B416 (1998) 36. 

[9] S. Ando and D.-P. Min, Phys. Lett. B 417 (1998) 177. 
[10] S. Ando, F. Myhrer and K. Kubodera, Phys. Rev. C 63 015203 (2001) 
[11] V. Bernard, T. Hemmert and U.-G. Meifiner, Nucl. Phys. A 686, 290 (2001) 

[12] H. W. Fearing, R. Lewis, N. Mobed and S. Scherer, Phys. Rev. D 56 (1997) 1783; Nucl. Phys. A 631 
(1998) 735c; V. Bernard, H. W. Fearing, T. R. Hemmert and U.-G. Meifiner, Nucl. Phys. A 635 (1198) 
121; A 642 (1998) 563. 

[13] D. D. Bakalov, M. P. Faifman, L. I. Ponomarev, and S. I. Vinitsky, Nucl. Phys. A 384 (1982) 302. 
[14] S. Ando, F. Myhrer and K. Kubodera, USC report USC(NT)-01-06, |nucl-th/0109"068| (2001) 
[15] S. Weinberg, Phys. Rev. Lett. 4, 575 (1960). 

[16] D. V. Bailin et al, PSI proposal R-97-05.1 (Dec. 1996); P. Kammel, priv. comm. and talk at the Intern. 

Riken Conference, MCF01, Shimoda, Japan (April 2001) to appear in Hyperfine Interactions 
[17] D. Armstrong et al, TRIUMF proposal E766 (June 1997). 



